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Assembly and Application of a Shock Data Analysis System 

1. Introduction 

Measured data from shipboard shock environments usually consist of 
analogue magnetic-tape recordings of acceleration and strain versus time. 
Usually, little analysis is required of strain recordings save spectral analy¬ 
sis and perhaps sums and differences. From the accelerometer recordings, 
however, it is usually necessary to extract velocities and displacements, and 
often shock spectra, as well. Shock data analysis facilities tend to run to 
two extremes of size and capacity. At one end is the large digital processing 
facility where hundreds of records (from an at-sea ship shock test, for 
example) are processed in batch, usually after the test is completed. At the 
under end is the small, often analogue, facility which can do limited analysis 
of a small number of records, although perhaps very slowly. The latter may be 
used to advantage for analyzing a few selected critical records in the field, 
or in the shock simulation lab. 

A need exists for a capability in the middle ground - a system capable of 
doing analysis on a modest number of records in a reasonable time, which 
accepts analogue magnetic tapes as input and produces report-ready graphs and 
tables as output, and which is field-transportable if necessary. 

This report describes such a system which has been assembled with 
standard, off-the-shelf instruments which are marketed nationwide and easily 
available. It is not in any way optimized for the specific task of shock data 
analysis, hence may not be competitive in speed, cost, etc. to a system which 
has been optimized in hardware and/or software to this specific application. 

It furnishes an example of what can be done along these lines with virtually 
no expenditure of R & D effort, and is an example in that any of its component 
instruments could be replaced by a different supplier's product with little 
change in operation or performance. 

2. Background 

2.1. Shock motions may be described in terms of acceleration -, 
velocity -, or displacement-time histories, and in all three parameters may 
vary greatly in waveform complexity. Simpler shock motions may be closely 
defined by a few waveform parameters, and these may also provide a reliable 
basis for comparison of shock motions. Shipboard shock, however, is extremely 
complex. Motions are typified by waveforms of high frequency range, high 
dynamic range and long duration which cannot be defined adequately by a few 
parameters. The additional descriptions found most useful for a motion of 
this type are the shock spectrum, which may be interpreted as a measure of the 
effects of a shock motion, and the design shock spectrum, which provides the 
data needed to design a structure to withstand it. Each provides a basis for 
comparing shock severities, and applies equally well to shock motions of all 
degrees of complexity of waveform. The general nature of shock motions, shock 
spectra and design shock spectra are discussed in Refs. (1) and (2). 

2.2. While a shock motion may be identified equally well by its acceler- 
Manuscript submitted May 29, 1981. 
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ation, velocity or displacement waveform, one of them will usually be most 
convenient for a given application. However, in most situations only acceler¬ 
ation can be measured properly. This is almost always the situation in the 
case of shipboard shock, and the bulk of data in this field will consist of 
records of acceleration and strain versus time, with perhaps a small fraction 
being records of velocity versus time. Thus, an analysis system for shipboard 
shock motion records must be able to produce time histories of acceleration, 
velocity and displacement, as well as shock spectra and design shock spectra, 
given input signals representing acceleration - or velocity-time histories, 
and preferably from displacement-time histories as well. Since shipboard 
shock motions are among the most demanding in terms of functional 
requirements, an analysis system capable of handling them will also be able to 
handle shock motions of most other types. 

3. Analysis System Characteristics 

3.1. The principal application of the system is the analysis of shock 
motions measured on laboratory shock machines. This includes those which 
generate simple-pulse shock, but most prominent are those which generate 
shipboard shock type environments: the Lightweight Shock Machine (LWSM), the 
Mediumweight Shock Machine (MWSM), the Floating Shock Platform (FSP) and the 
Large Floating Shock Platform (LFSP) (Ref. 3). Tests on these machines will 
be lightly instrumented compared to at-sea tests of ships and submarines, or 
some of the more specialized test vehicles. The typical instrunentation suite 
may be a few channels on the simpler-pulse machines and the LWSM, up to a 
dozen on the MWSM and FSP and a few dozen on the LFSP, while on some occasions 
several times as many channels may be required. Normally, about one-third of 
the total channels will fall in the category of equipment input motions, while 
most of the remainder will be equipment response motions and structural 
strains. 


3.2. The function of the analysis system is to accept analogue input 
signals, perform whatever processing is required, and produce output represen¬ 
tations of the desired quantities, scaled in the appropriate engineering 
units, in graphical or tabular form, and in video or hard-copy display. All 
signals require A-D conversion and storage, scaling and display. All motion 
signals require in addition computation and display of acceleration, velocity 
and displacement, each properly scaled, which requires the functions of 
integration and possibly differentiation. Input motion signals, and some 
response motion signals, also require the computation of maxiraax and residual 
shock spectra. This last computation is by far the most complex and time- 
consuming, but is needed for only a few channels. For a suite of a dozen 
channels of assorted information, the system will perform a complete analysis 
in an hour or two. 


Shock spectrum analyzers are available commercially, and are typically 
compact units and very fast. As stand-alone systems, they are lacking in that 
they cannot do a complete analysis (integration and scaling, etc.) and have no 
display. As components for calculating shock spectra they are wanting in that 
they cannot compute undamped spectra and their analysis frequencies are fixed. 
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3.3. Since Che analysis system is an assembly of component units designed 
for use as independent, general purpose instruments, it has a great degree of 
flexibility. With different switch settings, it would function equally well 
for analyzing earthquake motions or blast loadings. A brief description of 
each functional block of the system and its capabilities follows. 

3.3.1. Magnetic-Tape Recorder 

The primitive form of shock data records is usually FM 
magnetic tape. The recorder used with the analysis system has 14- 
track record/reproduce capability, plus edge tracks for voice annota¬ 
tion. The normal configuration has two housekeeping channels, one of 
them DR for a tape-speed servo reference signal and the other an FM 
channel devoted to time-code in IRIG B format, and 12 FM data 
channels. All FM channels are IRIG Intermediate Band I. Normal 
recording speed is 30 ips, giving a nominal bandwidth of 0 - 10 
KHz. The recorder can be operated at tapes speeds from 15/16 to 120 
ips, allowing the traditional analysis tricks of data speed-up and 
slow-down; record density can be changed by rearranging plug-in 
components on the electronics cards. Output level is normally set 
for lv_ rms for ±40 7 , peak-to-peak (sinusoidal) modulation. This 
particular unit was selected because it performs a set of canned 
calibration and set-up routines under control of a built-in micropro¬ 
cessor, and because it is a portable unit. Both features are 
convenient for field operations. 

3.3.2. Auxiliary Equipment 

The principal auxiliary equipment units are a time-code gener¬ 
ator/reader, a direct-writing oscillograph, and a storage CRO. In 
generator mode, the time-code unit produces an IRIG B code for 
recording in the magnetic tape. In reader mode, it monitors the 
recorded time-code and produces a pulse when a chosen time value 
occurs. This time can be selected to the millisecond and provides a 
start signal to the transient-storage unit (discussed below) and 
facilitates synchronization of the recorded data signals. The other 
units, the CRO and oscillograph, provide quick-look disolay of the 
recorded data signals and time-code. The CRO is a two-channel, 
single time-base type similar to many currently available. The 
oscillograph produces strip charts of up to 12 channels with a band¬ 
width of 0 - 5 KHz at paper speeds from 0.1 to 160 ips. 

3.3.3. Data Filter 

The data signals output from the tape recorder are passed 
through an 8-channel data filter. Each channel has a low-pass 
characteristic with cutoff frequency selectable from 1 Hz to 15 KHz, 
and is direct-coupled. The unique features of this unit are that the 
filter characteristic is optimized for linear phase-shift within 
pass-band, minimizing envelope distortion, and that phase-shift 
between channels is closely controlled at no more than 2° between any 
pair. The filter characteristic also features a cutoff rate of 43 
db/octave, so that aliasing is negligible for sampling frequencies of 
4.2 times the cutoff frequency or higher. Each channel has unity 
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gain, and can accept inputs of up to 10 _V . As used in the analysis 
system, the cutoff frequency is normally set to IKHz. 

3.3.4. Transient Storage Unit 

The transient-storage unit contains seven independent AD and 
digital storage channels with a common tine-base. Each channel can 
be operated in an inhibit mode which preserves the existing stored 
signal. The unit can be operated in a self-triggered mode, where the 
signal in any channel can be used to initiate the capture of the 
signals in all channels which are not inhibited. As used in the 
analysis system, the capture-initiate signal is derived from the 
time-code reader to preserve commonality of time-base between succes¬ 
sive batches of signals. Upon initiation, each channel captures 4096 
10-bit samples of its signal, for a dynamic range of 60 db. As used 
in the analysis system, the channel full-scale range is ± 2 V_ ; it 
can be varied from ±0.1 to ±5V_. The sampling period used is 100 
microsec., for a total record length of 400 millisec.; the sampling 
period can be varied from 5 microsec. to 2.5 sec. Each channel 
provides outputs of the digitized "staircase" and smoothed analogue 
versions of the stored signal waveform. As used in the analysis 
system, the signal stored in any selected channel is read out on the 
IEEE-488 bus in bit-parallel, character-serial ASCII code. This 
particular unit was selected because of its 7-channel capacity and 
the fact that each A-D converter and storage channel is functionally 
independent from the others. 

3.3.5. Processor and Display Unit 

The processor and display unit is a large-screen desk-top 
computer and peripherals. The computer has 64K-byte memory and a 
resident BASIC interpreter, leaving some 53K-bytes of usable 
memory. Additional ROM packs provide signal processing routines, 
including integration and differentiation, and text-editing 
capability. The principal display is a built-in large screen CRT. 

The peripherals include a hard-copy unit which gives a permanent 
record of the CRT display with reduced size, and a digital X-Y 
plotter which gives high quality graphs up to 11 x 17 in. Other 
peripherals are a high capacity disc-file manager, which allows a 
large library of programs to be resident on floppy disc, and an 
additional tape-cassette drive-unit, allowing more flexibility in 
data transfer. All of these components are produced by the same 
manufacturer, and are configured and interfaced to form a compatible 
system. This (processor and display) unit was selected because it 
possesses the requisite computational power and has an excellent 
repertoire of graphics routines. Other features are good man-machine 
interface, BASIC interpreter and large-screen display. These seem to 
assist the Intended user profile, engineering/technical personnel 
without extensive familarity with computing systems, in learning to 
operate this analysis system with a rainiraun of instruction. The 
existence of a comprehensive assortment of proven peripherals was 
also a factor of consideration. Another feature of the computer is 
that it can also be operated as a high-speed (9600 baud) intelligent 
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terminal, should operation with a larger host computer be 
desirable. The last peripheral in the system is a high-speed 
printing terminal from another manufacturer. When used with the 
analysis-.system, it functions as a moderately fast line-printer. 

3.4. The functional block diagram of the analysis system is in shown in 
Fig. 1. 

4. Computation of Shock Spectra 

4.1. Recalling that the shock spectrum of a shock motion may be 
interpreted as the maximum relative displacements of a set of massless single- 
degree-of-freedom systems whose bases are forced to follow the motion, it is 
apparent that shock spectra may be derived for any degree of damping. For the 
purpose of comparing shock motion environments and specifying design shock 
spectra, the undamped version is preferred, since the damped spectra are less 
sensitive to variations in the input motion. Shock spectra can be computed 
from the differential equation of motion, or from the corresponding integral 
equation. Either form can be expressed to involve the input motion in terms 
of acceleration, displacement or velocity. The procedure used here is that of 
numerical solution of the differential equation relating relative displacement 
to input accleration. The computational procedure requires the input to the 
program as the acceleration-time history of the shock motion. The 
computational procedure uses the algorithm described in Ref. (4). 

The equation of motion for an undamped oscillator of frequency o is 
(1) X + ui l X - -Z 


with the solution for the relative displacement X 

• . ** 

(2) X(t) = X Q cos wt +-^2. s in «t - — /£z(x) sin w(t - t) dx. 

The numerical integration equation derived by G. O'Hara (Ref. 4) is 


(3) 

X n+ ^w = X n u> cos wh + X n sin wh 


e 

n( 1 - cos wh) ri t /I + cos wh sin u>h\ 

wTi " s n-l ( -2SK-ST r F _) 


with a parabolic approximation for the digitized input function 
here S n - K+l - K 


Z (t); 


and S l n _i 
(note that S^ n _^ 
so that 


S n - Vl " Z n+1 - 2Z n 4- Z n _ x 
is the second difference, not a square) 


(4) 


Z(t) 
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In these equations, h is the time interval of integration. 

By eliminating X n in (3), R. Bort's implementation of equation 3 (Ref. 5) 
is 

(5) X(t + h) = 2 X(t) cos uh - X(t-h) + h [Z (t + h/2) - Z(t - h/2)] 

Equation (5) has been applied in this program because it consists of only 
three additions. As implemented here, the term ^n brackets 
(Z (t + h/2) - Z (t - h/2)] is identified as hZ (t), a sample point from the 
digitized input acceleration record. Before processing input data the 
2 cos uh term is determined and stored in array C. All scaling factors are 
deferred to the end of the calculation for each analysis frequency. 

In order to cut down calculation time, the fact that lower analysis 
frequencies require fewer periods to reach peak response than higher 
frequencies is employed. The number of periods per analysis frequency 
considered in the processing routine is given by 

(,;) 0 * nfi-F 

where F * w/2ir . 

The second procedure to speed up the program selects data points with respect 
to cutoff frequency and analysis frequency. The calculation time has been 
reduced this way from about 50 min. to 8 min. while producing only a minor 
error (see Fig. 2A,B). 

5. System Operator's Manual 

5.1.1. It is assumed that a time record of the shock motion - 
including a calibration step - is available on an analog tape and a 
time code is recorded on a separate track of the tape. Playing back 
the record (or calibration step) and the time code on the Visicorder, 
determine the millisecond preceding the starting time of the 
record. Set the time code generator to open the time gate and 
trigger the transient recorder. The transient recorder accomplishes 
the A/D - conversion and stores the digitized time record in its 
memory. An interface bus allows the computer to have access to the 
seven channels (memories) of the transient recorder. Further data 
processing is done under the control of the program listed in the 
Appendix. 

5.1.2. This program is written in BASIC and runs on the desktop 
computer Tektronix 4054. The three parts of the program are stored 
on the disc in three files (a back up tape cartridge contains the 
program in binary code). Part I reads the calibration step and the 
time record from the memory of the transient recorder, and calculates 
shock spectra. Part II contains all output procedures. Part III is 
basically the same as Part I except that it reads a calibrated 
digital time record from the external tape unit rather than from the 
transient recorder. Both Part I and Part III use the output 
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(a) Maximax shock spectrum 
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(b) Residual shock spectrum 

introduced by data decimation. Maximax (a) and residual (b) shock spectrum computed from the 
complete (4096) and decimated data sets for a half-sine acceleration pulse. 








procedures in Part II. 

5.2 Starting Up 

The computer and at least half of the units connected with the 
computer by interface bus must be turned on. To run the program, set the 
clock on the disc drive, mount the unit, load the program and enter "RUN". 
Example: Enter: CALL "SETTIM", "20 - MAR -81 □ 8:30" 

Press: RETURN 

Enter: CALL "MOUNT ",t , A$ 

Press: RETURN 

Enter: OLD " (? SPRCM/PART I" 

Press: RETURN 
Enter: RUN 
Press: RETURN 

Message on the Screen: 

SPRGM/PART I is loaded (default conditions)! 

The corresponding procedure for using binary tape is either: 

Enter: FIND 2 
Press: RETURN 
Enter: CALL "BOLD" 

Press: RETURN 
Enter: RUN 
Press: RETURN 
or simply press: AUTOLOAD 

While running the program, the computer indicates a need of communication 
with the user by a "beep" tone and a message on the screen. 

The symbol Q indicates one space. 

5.3 Default Conditions 

Default conditions are listed in line 290 to 500 in the Appendix. All 
default parameters may be changed at any time by entering the desired values 
from the keyboard (e.g. entering F2 ■ 500 changes the upper analysis frequency 
to 500 Hz). The current value of every parameter can be checked by entering 
the parameter and pressing RETURN: the current value of the parameter appears 
on the screen. 

The most important parameters defined by default conditions are explained 
below: 

FI, F2: Lower/upper analysis frequency: can be any 

positive value between 1 and 1000. 

C0: Number of analysis frequencies between FI and 

F2: depends on memory available (Part I up to 
about 110, Part III about 160). C0 * 1 is 
possible also. 

F8: Cut off frequency of low pass filter (default 

1000 Hz). 

Nl, N2: Lower/upper value of y-axis in log display 
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N6, N7: Lower/upper value of x-axis in log display 

D: Primary address for output devices: 

32 screen 

3 ext. tape 

4 plotter 
40 printer 

P: selects grid in log display: 

0<P<10 sets grid line at P 
P>10 no grid 
(P must not be zero!) 

0: indicates the part of the program currently 

loaded in the memory 

1. Part I 

2. Part II 

3. Part III 

(Do not change parameter 0 from the keyboard). 
M9: Calibration matrix. 

A2: Offset-correction matrix. 

5.4. User Defined Keys 


The first row of user defined keys (No. 11 through 15) are used with 
program Part I or III. All other user defined keys (No. 1 through 10 and 16 
through 20) belong to program Part II. User Defined Key No. 4 and _5_ allow 
switching from program Part I or III to Part II or vice versa. Pressing a 
user defined key which belongs to a program part not currently loaded in the 
memory causes an error message to appear on the screen. 

The following describes each user defined key and the system's 

reaction. 


5.4.1. First Row (No. 11 through 15) 

User Defined Key No 11 triggers the calibration routine in program 
Part I. After the time-record has been defined as a displacement, 
velocity or acceleration record there are two ways to fill the 
calibration matrix M9 with the calibration factors for each channel. 

First the calibration factor may be entered from the keyboard in input 
units per 1 Volt, or secondly, the calibration step may be read from each 
channel of the transient recorder. In the latter case the calibration 
level has to be entered in appropriate units before the computer reads 
the content of each channel and displays it on the screen. The user may 
define the calibration level by means of setting the crosshair first to 
the upper and then to the lower calibration level. The calibration 
factor (i.e. input unit/digit) for this particular channel will be 
assigned to an element of calibration matrix M9. The program starts this 
procedure all over again until the calibration steps of all 7 channels 
have been read. If the program encounters a channel set to ground or 
calibration mode an appropriate message appears on the screen and no 
calibration factor will be defined for this particular channel. The 
program proceeds with the next channel. 

User Defined Key No. 12 reads the time record from a defined 
channel of the transient recorder and displays a plot of acceleration. 
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velocity and displacement's time on the screen. Also each maximum and 
minimum value is printed. When running program Part I, the calibration 
routine (key 11) must be completed before using Key No. 12 . With program 
part III, User Defined Key No. 12 causes the computer to read the time 
record (including other important information like calibration factor, 
time step, etc.) from the external tape. 

The read routine has to be completed for each record before 
calculating shock data by pushing User Defined Key No. 13. The computer 
then calculates shock data of the time record that has been input by 
means of Key 12. 

Using the mathematical algorithm described above, the program 
computes: 

1. Analysis Frequency 

2. Pseudo velocity - Maxiraax-Shock Spectrum 

3. Rel. Velocity - Residual-Shock Spectrum (or 
Fourier Transform) 

4. Rel. Displacements 

5. Abs. Acceleration 

6. Time when maximum occurs 

7. Phase 

While calculating, the computer displays on the screen the serial 
number of analysis frequency (from 1 at the lowest to Zb at the 
highest), the analysis frequency, the pseudo-velocity, and the 
number of data points used for the calculation of shock spectra for 
each analysis frequency. The computing time is also given at the 
end of the procedure. 

The number of analysis frequencies (1 to 160 ), as well as any 
interesting range of analysis frequencies between 1 and 1000 Hz, can 
be chosen. (Note that F2 has to be larger than Fj). 

The time record and all important information, such as 
calibration factor, time step, etc. can be stored on a defined file 
of the external tape by pushing User Defined Key No. 14. Additional 
analysis can be done on this time record by means of part III of the 
program (refer to User Defined Key No. 12). An offset correction of 
the time record can be made with the procedure of User Defined Key 
No. 15. The first integral of the time record is displayed on the 
screen and the user may define the zero line by means of the cross 
hair. After pushing any key the computer calculates an offset 
correction for this particular channel and considers that at all 
later operations. This procedure may be repeated several times. 

* The title or description of the time record is stored in string W$ and 
appears on the top of the display. It may be changed anytime by redefining W$ 

^Depends on memory available for the calculation (program part III allows up 
to 160 analysis frequencies). 
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However, the user has to realize that this procedure requires some actual 
value to be known. 

5.4.2. Second Row (No. 1 through 5) 

The results of the shock spectra computation ( Key No. 13) can 
be saved by pushing User Defined Key No. 1 . This causes the 
computer to write all data on a defined file of the external tape. 
These shock data can be read back with User Defined Key No. 2 . 

A numerical output of all analysis data can be retrieved by 
pushing User Defined Key No. 3 . As for all other output procedures 
the primary address of the selected output device (screen, printer, 
or plotter) has to be specified here. 

User Defined Key No. 4 makes the computer load part II of the 
program while keeping all values of all parameters and variables in 
its memory. 

5.4.3. Third Row (No. 16 through 20) 

The third row of user defined keys produces a linear plot of 
the shock sgectra. The amplitude scale is automatically adjusted to 
the output. Any range of analysis frequencies can be selected. 

User Defined Key No. 16 displays a plot of pseudo-velocity 
vs. analysis frequency (maximax) 

User Defined Key No. 17 displays a plot of rel. velocity 
vs. analysis frequency (Residual) 

User Defined Key No. 18 displays a plot of acceleration vs. 
analysis frequency. 

User Defined Key No. 19 displays a plot of rel. 
displacements vs. analysis frequency. 

User Defined Key No. 20 displays a plot of phase (0 to 180 
degree) vs. analysis frequency. 

5.4.4. Fourth Row (No. 6 through 10) 

A logarithmic plot of analysis data can be retrieved by 
pushing any of the user defined keys in the fourth row. User 
Defined Key No. 6 displays the Maximax Shock Spectrum, No. 7 the 
Residual Shock Spectrum, and No. 8 both of them. User Defined Key 
No. 9 allows some analysis work on the screen in determining modal 
frequencies and appropriate shock spectrum values (the Design Shock 
Spectrum) by placing the two crosshairs at points of interest. The 
thumbwheels on the keyboard control the movement of the crosshairs 
When a keyboard key is pressed, the location of the crosshair is 
marked and the appropriate analysis frequency and shock spectrum 
value is printed beside the graph-plot. When the key "0" is 


* 


except phase is always displayed between 0 and 180 degrees. 
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pressed, the crosshair disappears from the screen. (Note: this 
procedure works with the screen only!) 


User Defined Key No. 10 gives a plot of a number of different 
shock spectra stored on the external tape. Fach spectrun is drawn 
with a different dark line and the description of the spectrum is 
printed wherever the crosshairs have been placed. 

6. Examples of Analysis System Output. 

6.1. Figures 3 and 4 show the various output formats which can be 
selected. These were made using the hard copy unit, which produces a 
reduced-size version of the CRT display. Higher quality graphs can be 
obtained from the plotter. 

6.2. For Fig. 3, the input record was a half-sine acceleration pulse. 

The original, integral and double integral are shown in Fig. 3A as scaled 
acceleration, velocity and displacement. The maximax shock spectrum is 
shown in linear plot, without grid-lines, in Fig. 3B, and in log-log 
plot, with grid-lines, in Fig. 3C. The maximax and residual are shown in 
log-log plot, with factor-of-ten grid-lines, in Fig. 3D. Finally, a 
tabular listing of all computed data is shown in Fig. 3E. 

6.3. The input record for Fig. 4 was an acceleration-time history 
recorded on a resilient structure aboard the Floating Shock Platform with 
shot standoff of 20 ft. The acceleration, velocity and displacement 
waveforms are shown in Fig. 4A, utilizing the offset correction routine 
described above: i.e., requiring the average velocity near the 400 
millisec. mark to be zero. The maximax and residual spectra are given in 
Fig. 4B, and again in Fig. 4C, which demonstrates the extraction of 
design shock spectral parameters. Modal frequencies and shock spectrum 
values are lined up with the CRT crosshairs, and their values tabulated. 
The values read are also identified on the graph. The tabular listing of 
computed data is shown in Fig. 4D. 
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11 m* record (a) Acceleration, velocity and displacement time histories 

Fig, 3 — Output for analysis of a half-sine acceleration pulse 
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Maximax shock spectrum, linear plot, without grid 
— Output for analysis of a half-sine acceleration pulse 
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(d) Maximax and residual shock spectra, log-plot, with factor-of-ten grid 
Fig. 3 — Output for analysis of a half-sine acceleration pulse 
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(e) Tabular listing of all computed data 
Fig. 3 — Output for analysis of a half-sine acceleration pulse 
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(a) Acceleration, velocity, and displacement time histories 
Fig. 4 - Output for analysis of FSP acceleration record 











(Zero Doaofrg) Noiurol Frequency Hz 

(b) Maximax and residual shock spectra 
Fig. 4 — Output for analysis of FSP acceleration record 
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(d) Tabular listing of all computed data 
Fig. 4 — Output for analysis of FSP acceleration record 
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Appendix 1 

List of Program Variables 


1 A Arroy for rime record OITKN0) 

2:A2 Zero offset - Matrix 010(7,2) 

3 B Array 010(21 

4 B3 Voriadle used in sudrouiine TIOE 

5 84 vorioole used m subroutine TIOE 

6 C Arroy OIIKC0) 

7;C0 Actual numoer of analysis frequencies (defoult 90) 

8:CI Element of F5 (firsi analysis frequency per frequency range) 

9:C2 Element of F5 (losi analysis frequency oer frequency ronge) 

10.C3 Max possible analysis frequencies (depends on memory availoble) 

110 1/0 device oddness (defouli 32) 

12:0! -i1e number on ext. tape drive to store data 

1303 File number on ext. tope drive to read data from 

14.E Instantaneous factor for sample step 

15 El "actor for sample siep within frequency ronge 1U-10HZ) 

16 E2 Factor for somple siep within frequency range il(10-100Hz> 

17E3 Factor for somple step within frequency ronge 111(!00-1000Hz! 

1 SF Instantaneous value of analysis frequency 

19F1 Lowest analysis frequency (default 1) 

20:F2 Highest analysis frequency 1 default 1000) 

21 :F3 -owe'" analysis frequency wiihin one of the 3 frequency ranges 

22:F4 Upper analysis frequency within one of the 3 frequency ranges 

23:F5 Frequency array 01HICB) 

24 : F6 Lower frequency of lin. graph (input) 

25:F7 Upper frequency of lin. graph (input) 

26:F8 Cutoff frequency of lowpass filter 

27 G Instantaneous value of accelaration ((2*pi?F)f2*X) 

28 G5 Acceleration array OI(1(C0) 

29 H Time increment of digitized time record 

30 HI Exponent 

3" H2 Time increment for computation 

32 I Counter 

33 It Courier to find Cl 

34-12 Indicates frequency range (I to III) 

3513 Sample number when V2 occurs 
36.14 Sample number when VI occurs 
37 K Target or-ay for graph disp lay 

38'M Factor of F8 to calculate somple frequency withm frequency ronae I Id 
e'aul’ 2) 

39 *2 Factor of F8 to calculate somple frequency withm frequency range II ( 
defaulr 2) 

40 *3 Factor of F8 to calculate sample frequency withm frequency ronge II l 
idef ault l0) 

41 L Calibration level in input units 
42.F Number of data module 

43 HI Maximum of arroy N for autoscole 

44 F9 Calibration Hoirix 0IM(7,2) 

45 N0 Number of data points on record (default 4096) 

46 N Njmber of data points considered during calculation 

47 ni _ower vertical level of log graoh display 

48 N2 Upper vert.col level of log graph display 

49 N4 Number of data points at the beginning of rne residual response 
50-N5 Number of data points after trunctotion 

5' N6 Lower frequency of log graph display 
52 N7 uoper frequency of log graph display 

53.C I no 1 cates *he port of the program currently loaded M or 2) 

54 P Selects g-•d on log. display (P>!0 no grid, defoult 10) 
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55.PI Error indicator if >0 902 

5S Q Number of periods considered per analysis frequency 

57.R Selects (lAXIOAX- and/or RESIDUAL SHOCK SPECTRuflt Keys S to 9l 

58 R1 Selects div. linear spectra (user defined Keys IS to 20) 

59 R2 Ind cates displacement(!).velocity(2) or occelerotion(3> record 

60 R3 Equals X3 before N4 (initial I 
61.R4 Equels X3 after NA (residual! 

62 R5 Flag for offset-procedure 

63 R6 Flag for deleting colifcrotion routine 

64 R 8 Flag m col Brotion routine 

65 S CaliDrat ion factor 

66 SI Calibration factor after Oifferern iotion and/or integration 

67 ■ 62 Calibration factor after E-procedure 

63 T 1 nstontaneous value of time of maximum (T3 or T 4 ) 

69-Ti Beginning of time record 

70 T3 Time of max. response before N4 'initial* 

7' 74 Time o‘ max. response after N4 (residual) 

72 T5 Time array (contains max values of either T3 or T41 Oin(C0i 

73 U Time difference m aecimal value 

'•’4 V Instontoneous value of pseudveloc i ty ((2*pi*F)*X) 

75:V t Maximum of orroy A 
78 V2 Mini mum of array A 

77 V5 Velocity array OIM(C0) 

78 : if Factor 

7g u' Number of spectra pro plot (Multi Display node) and flag for ws-qispla 

y 

30 W2 Array. Old W2(WI) 

81 X X-coordmate of pointer 

82 X Instantaneous value of displacement (maximum of R3 or R4) 

S3 X- Deflection at -0.5*H2 

84 X2 Oeflection at +0.5SH2 

65:x3 Max. response of analysis frequency F 

86:X5 Displacement array 

87 v '-coordinate of pointer 

88 y Instanraneous value of of max. response velocity after N4 
8S:Y5 Veloc.ty array OIMIC0! 

90 Z Instantaneous value of of phase 
Si:Z5 Phase array D[fl(C0) 

92:as String for graph display 

93:3* System time when process-procedure is completed 
94:C* T,me Siring 

95.0* Siring to address data module 
96 G* String far g^aph display 
97.H* String far graph display 

98 l* String -0 checK cond.’ion of data module 

99 MS String for graph display 

1 00 ns String for grooh display 

101 PS Spring far graph display 

!02 Qt Spring t f '0" slops anolysis procedure) 

03:R* "NRl - CCOE 5842" 

104 TS Contain* system time 
05:Vs Contain* text for vert level 
' 06 Vi Title to identify record 
1 07:x* Segment string 
'08 ys Segment strmg 
109 Z* Segment string 
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Appendix 2 
Program Listing 

1 REM //////////////////////// 9SPRGM/PART I \\\\\\\\\\\\\\\\\\\\\\\\\\\ 

2 INIT 

3 GO TO 260 

4 IF 0=1 THEN 92 

5 GO TO 7450 

8 IF 0=1 THEN 92 

9 GO TO 7530 

12 IF 0=1 THEN 92 

13 GO TO 8980 
16 GO TO 190 
20 GO TO 100 

24 IF 0=1 THEN 92 

25 R*1 

26 GO TO 7690 

28 IF 0=1 THEN 92 

27 R =2 

30 GO TO 7690 

32 IF 0=1 THEN 92 

33 R=3 

34 GO TO 7690 

36 IF 0=1 THEN 92 

37 GO TO 9510 

40 IF 0=1 THEN 92 

41 GO TO 9780 

44 IF 0=2 THEN 90 

45 CHARSIZE 3 

46 IF R6=0 THEN 94 

47 IF R6=l THEN 4200 

48 IF 0=2 THEN 90 

49 GO TO 530 

50 DELETE A 

52 IF 0=2 THEN 90 

53 GO TO 2510 

54 DELETE A 

56 IF 0=2 THEN 90 

57 GO TO 5250 

60 IF 0=2 THEN 90 

61 R5=l 

62 CHARSIZE 3 

63 GO TO 2410 

64 IF 0=1 THEN 92 

65 R1=1 

66 GO TO 6010 

68 IF 0=1 THEN 92 

69 R1=4 

70 GO TO 6010 

72 IF 0=1 THEN 92 

73 R1=2 

74 GO TO 6010 

76 IF 0=1 THEN 92 

77 R1=3 

78 GO TO 6010 

80 IF 0=1 THEN 92 

81 R1=5 

82 GO TO 6010 

90 PRINT "You qot to PRESS KEY 5 first!G" 

91 RETURN 

92 PRINT "You qot to PRESS KEY 4 firstlG" 

93 RETURN 

94 PR I "LJJTHE CALIBRATION PROCEDURE HAS BEEN DELETED IN ORDER TO SAUE"; 
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95 PRINT " MEMORY. IF YOU WANT TO REDEFINE THE CALIBRATION MATRIX M9 

96 PRINT "RELOAD THE PROGRAM (PRESS USER DEFINED KEY 4 AND 5 IN A ROW" 

97 PRINT "OR PRESS >AUTOLOAD<) 1 GG" 

99 RETURN 

100 REM *************** LOAD ROUTINE *********************************** 
110 IF 0=1 THEN 160 

120 DELETE F5,X5,05,G5,T5,Y5.Z5,C,K , A 
130 R6=1 
140 0=1 

150 APPEND "8SPRGM/PARTI" ; 1 

160 PRINT “LGSPRGM/PART I IS LOADED!” 

170 DELETE 5380,20000 
130 RETURN 
190 IF 0=2 THEN 240 
200 0=2 

210 DELETE 270.20000 

220 APPEND "8SPRGM/PARTII";260 

230 RENUMBER 6000,10,260 

240 PRINT "LGSPRGM/PART II IS LOADED!” 

250 RETURN 

260 REM *************** DEFAULT CONDITION ****************************** 
270 Wt= "TEST" 

280 0=1 
290 W 1 = 0 
300 R5=0 
310 R6=1 
320 D=32 
330 F1=1 
340 F2=1000 
350 F8=1000 
360 N1=0.01 
370 N2 = l00 
380 N6=1 
390 N7=1000 
400 C0=90 
410 N0=4096 
420 K1=2 
430 K 2 = 2 
440 K3 =10 
450 P=10 

460 DIM M9(7,2),C(C0),A2<7,2),R$(15),Tt<18),Bt(18),Ct(18) 

470 ^IM Dt<8>,Ht<20>,Vt(22>,Gt(7>,M*(7),Nt(8>,X$(22>,P*(35>,Lt(4>,Qt(1> 
480 R4="NRL - CODE 5842" 

490 M9-0 
500 A2= 0 

510 PRINT " GSPRGM/PART I IS LOADED (DEFAULT CONDITIONS)!” 

520 RETURN 

530 REM ****************** READ RECORD ********************************* 
540 DELETE 4200,5240 
550 R6=Q 

560 PRINT 80,17:1.2,1.3 
570 CHARSI2E 3 

580 PRI "LTHE CURRENT TITLE IS '" ;Wt;”' J ENTER A NEW TITLE IF NECESSA”; 
590 PRINT "RY. OTHERWISE PRESS 'RETURN',JG" 

600 INPUT At 

610 IF At="" THEN 630 

620 Wt =AS 

630 DIM A < NO) 

640 CHARSI2E 3 
650 IF R5=0 THEN 660 
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660 W3*130 

670 IF D< >4 THEN 690 
680 W3=150 

690 PRINT "JENTER NUMBER OF DATA MODULE TO READ FROM,• G "; 

700 INPUT M 

710 IF R3=l THEN 740 

720 IF M9<M,2>=0 THEN 2390 

730 S=M9< M, 2 > 

740 PRINT 86!”P" 

730 GO TO M OF 760,780,800,820,840,860,880 
760 D«="I01B4096" 

770 GO TO 890 
780 D*="I02B4096" 

790 GO TO 890 
800 D*="I03B4096" 

810 GO TO 890 
820 D* a "I04B4096" 

830 GO TO 890 
840 D4= "105B4096" 

850 GO TO 890 
860 D<=”I06B4096” 

870 GO TO 890 
880 D*="I07B4096" 

890 PRINT e6;D* 

900 INPUT 86:A 
910 p*="D " 

920 PRINT 86: P* 

930 INPUT 86: P* 

940 A=A-A2(M,2) 

950 HS- “SA” 

960 PRINT 86:H* 

970 INPUT 86:H* 

980 HI=VAL(H*) 

990 IF HI>12 THEN 1180' 

1000 GO TO HI OF 1,1,1,1010,1030,1050,1070,1090,1110,1130,1150,1170 

1010 H=0.005 

1020 GO TO 1200 

1030 H=0.01 

1040 GO TO 1200 

1030 H=0.025 

1060 GO TO 1200 

1070 H=0 . 05 

1080 GO TO 1200 

1090 H=0.1 

1100 GO TO 1200 

1110 H=0.25 

1120 GO TO 1200 

1130 H-0.5 

1140 GO TO 1200 

1150 H=1 

1160 GO TO 1200 

1170 H=2.5 

1180 PRINT "ENTER TIME INCREMENT IN MS "; 

1190 INPUT H 

1200 IF R5=1 THEN 2430 

1210 CHARSIZE 2 

1220 CALL "TIME",T* 

1230 PRINT 8D: "LK";Rt 
1240 PRINT 8D:T$ 

1250 UINDOU 0.U3,0,100 
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1260 VIEWPORT 0,U3,0,100 
1270 IF W1>0 THEN 1350 
1280 PRINT 8D,17:2,3 
1290 CHARSIZE 4 
1300 MOVE 8D:W3/2,98 
1310 FOR 1=1 TO L£N<W*>/2 
1320 PRINT @D:"H"; 

1330 NEXT I 
1340 PRINT @D:W* 

1350 PRINT SD,17!1.2,1.8 

1360 CHARSIZE 2 

1370 T1=H#ABS(VAL<P4> > 

1380 GO TO R2 OF 1390,1520,1660 

1390 GOSUB 1840 

1400 S1=S 

1410 GOSUB 2020 

1420 PRINT @Di“*" 

1430 CALL "DIF3" , A > A 

1440 S1=S1/H 

1450 GOSUB 1900 

1460 GOSUB 2020 

1470 CALL "DIF3",A,A 

1480 S1=S1»10 A 3/<H*9.8044> 

1490 GOSUB 1960 
1500 GOSUB 2020 
1510 GO TO 1780 
1520 GOSUB 1900 
1530 S1=S 
1540 GOSUB 2020 
1550 PRINT @D : 11 *" 

1560 CALL "DIF3" , A, A 
1570 S1=S1*10 A 3/<9.8044*H) 

1580 GOSUB 1960 

1590 GOSUB 2020 

1600 CALL "INT ",A,A 

1610 CALL "INT ",A,A 

1620 S1=S1*1.OE-3»H a 2*9.8044 

1630 GOSUB 1840 

1640 GOSUB 2020 

1650 GO TO 1780 

1660 GOSUB 1960 

1670 S1=S 

1680 GOSUB 2020 

1690 PRINT @D; 

1700 CALL "INT",A,A 

1710 S1=S1*1.0E-3*H»9.8044 

1720 GOSUB 1900 

1730 GOSUB 2020 

1740 CALL "INT",A,A 

1750 SI=S1*H 

1760 GOSUB 1840 

1770 GOSUB 2020 

1780 VIEWPORT 0,W3,0,100 

1790 WINDOW 0,W3,0,100 

1800 MOVE @D:0,0 

1810 PRINT @D: "* oriq. ti*e record 
1820 DELETE A 
1830 RETURN 

1840 VIEWPORT 5,U3-30,4,32 
1850 V*=“DISPLACEMENT" 





1860 G*="D Max : " 

1870 M*="D Min 
1880 N*=“ [mm]" 

1390 RETURN 

1900 VIEWPORT 5,W3-30,36,64 
1910 V*="VELOCITY" 

1920 G*="V Max : " 

1930 M$="V Min : " 

1940 N$=" Im/sbcI" 

1950 RETURN 

1960 VIEWPORT 5,W3-30,68,96 
1970 V*="ACCELERATION" 

1980 G*="A Max :" 

1990 M*="A Min :" 

2000 N*=" [g]" 

2010 RETURN 

2020 CALL "MAX",A,VI,14 
2030 CALL "MIN",A,V2,13 
2040 WINDOW 0,NO,V2,VI 
2050 IF D=32 THEN 2110 
2060 MOVE 8D:0,0 
2070 FOR 1=1 TO NO 
2080 DRAW @D:I,A<I) 

2090 NEXT I 

2100 GO TO 2120 

2110 CALL “DISP " , A 

2120 AXIS @D:400,0 

2130 IF R5=l THEN 2340 

2140 FOR 1=0 TO NO STEP 2000 

2150 MOVE @D:1,0 

2160 PRINT @D:"HJ",I#H; 

2170 NEXT I 

2180 PRINT @D:"HHHHHHHHHHH[MSI" 

2190 MOVE @D:4200,V1 

2200 PRINT @D: USING 2240:"J",G«,V1*S1,N*j 
2210 MOVE @D:420 0 ,VI 
2220 PRINT @D:"J"; 

2230 PRINT @D: USING 2240i"J",M»,V2*S1,N$; 

2240 IMAGE A,7A,6D.1D,11A 
2250 MOVE @D:0,0,5*(V2+V1) 

2260 PRINT @D:"HH" ; 

2270 FOR 1=1 TO LEN<V*>/2 
2280 PRINT 8D: "K"; 

2290 NEXT I 

2300 FOR 1=1 TO LEN<V*> 

2310 X<=SEG(V$,1,1) 

2320 PRINT @D;Xi;"HJ"; 

2330 NEXT I 
2340 RETURN 

2350 PRINT @D:"JDATA MODULE ";M>" IS IN GROUND MODE MGG" 
2360 RETURN 

2370 PRINT "JDATA MODULE ";M;" IS IN CALIBRATION MODE MGG" 
2380 RETURN 

2390 PRINT "JREAD THE CALIBRATION STEP FIRST'S' 

2400 RETURN 

2410 PRINT "LOFFSET - CORRECTION!" 

2420 GO TO 630 
2430 VIEWPORT 0,110.10,90 
2440 CALL "INT",A,A 
2450 GOSUB 2020 
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2460 POINTER X,Y,Q* 

2470 A2(M,1)=M 

2480 A2<M,2)=Y/X+A2<M,2> 

2490 R5=0 
2S00 RETURN 

2510 REM a##**#*****#******## PROCESS a********************************* 
2520 CHARSIZE 3 
2530 CALL "TIME",B* 

2540 PRINT @D:"L"j B$;"J" 

2550 PRINT @32,26:2 
2560 IF FI<1 THEN 2620 
2570 IF FI>1000 THEN 2620 
2580 IF F2<1 THEN 2630 
2590 IF F2>1000 THEN 2630 
2600 IF F2< =F1 THEN 2630 
261« GO TO 2660 
2620 IF C0=1 THEN 2660 

2630 PRI “ANALYSIS FREQUENCY IS OUT OF RANGE ENTER NEW VALUES OF F1,F2G” 
2640 INPUT F1,F2 
2650 GO TO 2530 

2660 DELETE F5,X5,V5,G5,T5,Y5,Z5,C, K ,A 
2670 DIM A(NO) 

2680 C3=INT< <MEMORY/8-18)/8> 

2690 DIM F5<C0),X5<C0 >,V5<C0 >,G5(CO),T5<C0>,Y5(C0>,Z5(CO),C(C0) 

2700 £ 1=INT< 1000/<K1*F8*H)) 

2710 E2=INT< 1000/< K2#F8#H)) 

2720 E3=INT(1000/<K3#F8*H)> 

2730 F5< 1)=F1 

2740 IF C0=t THEN 2780 

2750 FOR 1=2 TO CO 

2760 F5<I)=<F2/F1> A (1/<C0-1>)*F5<I-1> 

2770 NEXT I 
2780 X5=0,01 
2790 V5=0,01 
2800 G5=«Q , 01 
2810 T5=0.01 
2820 Y5=0.01 
2830 Z5=0.0l 
2840 11=1 
2850 12=0 

2860 IF FI<10 THEN 2900 
2870 12=1 

2880 IF F1<100 THEN 2900 

2890 12*2 

2900 GO TO 3070 

2910 E=E1 

2920 GOSUB 3160 

2930 IF II>C0 THEN 3010 

2940 GO TO 3070 

2950 E*E2 

2960 GOSUB 3160 

2970 IF I1>C0 THEN 3010 

2980 GO TO 3070 

2990 E*E3 

3000 GOSUB 3160 

3010 DELETE A 

3020 CALL "TIME" ,T* 

3030 GOSUB 4030 
3040 PRINT 932,26:0 

3050 PRINT 9D:"_MAX. POSSIBLE ANALYSIS FREQUENCIES: ">C3 
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3060 RETURN 
3070 F3=F5(I1) 

3080 Cl=It 

3090 12=12+1 

3100 FOR 11=1 TO CO 

3110 IF F5<II)>10 A I2 THEN 3130 

3120 NEXT II 

3130 F4=F5(11 — 1> 

3140 GO TO 12 OF 2910,2950,2990 

3150 RETURN 

3160 PRINT 86: "P";D$ 

3170 INPUT 86:A 
3180 A=A-A2<M,2) 

3190 IF E<=1 THEN 3260 
3200 J=1 

3210 FOR 1=1 TO NO STEP E 
3220 A(J)=A<I) 

3230 J*J+1 
3240 NEXT I 
3250 GO TO 3270 
3260 E=1 
3270 N5=N0/E 
3280 H2=H*E 

3290 GO TO R2 OF 3300,3340,3370 

3300 CALL "DIF3" ,A,A 

3310 CALL "DIF3 " ,A , A 

3320 S2=S 

3330 GO TO 3380 

3340 CALL "DIF3",A,A 

3350 S2=S*H2 

3360 GO TO 3380 

3370 S2=S*10 A -3*H2 A 2»9.8044 

3380 C2=I1-1 

3390 P1=0.0062831853*F3*H2 
3400 FOR I=C1 TO C2 
3410 C<I)*2»C0S<P1) 

3420 IF Cl=C2 THEN 3440 
3430 P1=P1»(F4/F3)*(1/(C2-C1)) 

3440 NEXT I 

3450 IF P1< 0.902 THEN 3500 
3460 IF 143.566/H2>F2 THEN 3500 

3470 PRINT “LGUARNING: FREQUENCIES ARE TOO HIGH FOR YOUR TIME INCREMENT- 
3480 PRINT " ERRORS ARE OUTSIDE THE RANGE -7 TO *3 PERCENT FOR" 

3490 PRINT ” FREQUENCIES ABO^E ";143.566/H2;" HERTZ- 

3500 FOR I=C1 TO C2 
3510 F=F5(I) 

3520 IF F<=1 THEN 3550 
3530 N=1 000/(H2*5#LGT <F ) ) 

3540 IF N<N5 THEN 3560 

3550 N=N5 

3560 X1=0 

3570 X2=0 

3580 R3*0 

3590 R4=Q 

3600 IF F=0 THEN 3630 
3610 N4=N-1-500/< F *H2 > 

3620 IF N4>3 THEN 3640 
3630 N4*3 

3640 FOR J=1 TO N-l 
3650 X3=C (I)*X2-X1+A < J) 
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3660 IF J<N4 THEN 3710 

3670 IF ABS(X3>(ABS(R4) THEN 3740 

3680 R4=X3 

3690 T4=J 

3700 GO TO 3740 

3710 IF ABS (X3 > ( ABS(R3) THEN 3740 

3720 R3=X3 

3730 T3=J 

3740 XI=X2 

3750 X2=X3 

3760 NEXT J 

3770 IF ABS(R3)>ABS(R4) THEN 3800 

3780 R3=R4 

3790 T3=T4 

3800 X=S2*R3 

3810 U=0.00628318*F*X 

3820 G=6.288318»F*V/9.8044 

3830 T=Tl+(T3-0.5)*H2 

3840 Y=Q.00628318*F*S2*R4 

3850 Z=1.0E-3*(T1+(T4-0.5)*H2) *F + 0,25 

3860 Z=Z-INT(Z) 

3870 IF Y > 0 THEN 3890 
3880 Z=Z-0.5 

3890 IF Z<=0.5 THEN 3910 

3900 Z=Z-1 

3910 Z=360*Z 

3920 F5<I)=F 

3930 XS(I)=X 

3940 U5(I)=V 

3950 G5(I)=G 

3960 T5(I)=T 

3970 Y5(I)=Y 

3980 Z5(I)=Z 

3990 PRINT @Di USING 4000:I,F,V,N> 

4000 IMAGE 20T,3D,4X,4D.1D,4X,5D.2D,3X,6D 
4010 NEXT I 
4020 RETURN 

4030 REM ****** SUBROUTINE TIME ************************************** 
4040 DIM B(2) 

4050 C*=B* 

4060 FOR 1=1 TO 2 
4070 X*=SEG(CS.11,8) 

4080 Y*=SEG(X*,4,2) 

4090 Z$=SEG(X«,7,2) 

410 0 B(I)=UAL(X*) + (UAL(Y*)*60 +VAL < Z% ) >/360 0 
4110 C*=T* 

4120 NEXT I 
4130 U*B(2)-B(1) 

4140 B3=(U-INT(U))*60 
4150 B4=(B3-INT(B3>)*60 

4160 PRINT @D:"_CALCULATION TIME (hr :Min:sec): ■; 

4170 PRINT @D : USING 41 80 : INT (U> *MNT( B3 INT (B4) 

4180 IMAGE 2D,A,2D,A,2D 
4190 RETURN 

4200 REM *************** CALIBRATION *********************************** 
4210 DIM A (NO > 

4220 CHARSIZE 3 

4230 PRINT "LARE YOU DEALING WITH J A DISPLACEMENT RECORD (MM)''; 

4240 PRINT " (ENTER 1) A VELOCITY RECORD (M/SEC) (ENTER 2) OR A”; 

4250 PRINT - ACCELERATION RECORD (G) (ENTER 3)g "; 
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4260 INPUT R2 

4270 PRI " DO YOU WANT TO READ A CALIBRATION STEP (ENTER 1) OR"; 
4280 PRINT " TO ENTER CALIBRATION FACTOR (ENTER 2>G "; 

4290 INPUT R8 
4300 M9 = 0 

4310 IF R8=2 THEN 4960 
4320 FOR M=1 TO 7 

4330 PRI "LENTER CALIBRATION LEVEL FOR CHANNEL ";H;” IN INPUT UNITSG:"; 

4340 INPUT L 

4350 GOSUB 4750 

4360 IF S4=0 THEN 4390 

4370 IF S4=7 THEN 4420 

4380 GO TO 4450 

4390 GOSUB 2310 

4400 CALL "WAIT",2 

4410 GO TO 4730 

4420 GOSUB 2330 

4430 CALL "WAIT",2 

4440 GO TO 4730 

4450 PRINT @6;"P" 

4460 GO TO M OF 4470,4490,4510,4530,4550,4570,4590 
4470 D4= "101B4096" 

4480 GO TO 4600 
4490 D$="I02B4096“ 

4500 GO TO 4600 
4510 D*=“I03B4096“ 

4520 GO TO 4600 
4530 D$="I04B4096" 

4540 GO TO 4600 
4550 D*="IQ5B4096" 

4560 GO TO 4600 
4570 D*="I06B4096" 

4580 GO TO 4600 
4590 D$="I07B4096" 

4600 PRINT @6:D* 

4610 INPUT @6:A 

4620 WINDOW 0,5000,-512,512 

4630 VIEWPORT 0,130,0,100 

4640 PRINT "LCHANNEL ";M;" JJSET CROSSHAIR TO UPPER CALIBRATION LEVELG" 
4650 CALL "DISP",A 
4660 AXIS 

4670 POINTER XI,Y1,Q« 

4680 PRINT " A JJJ AND LOWER CALIBRATION LEVELG!" 

4690 POINTER X2,Y2,Q* 

4700 IF Y2>Y1 THEN 5230 

4710 M9<M,1)=M 

4720 M9(M,2)*L/(Y1-Y2) 

4730 NEXT M 
4740 RETURN 

4750 GO TO M OF 4760,4780,4800,4820,4840,4860,4880 
4760 L*="A01" 

4770 GO TO 4890 
4780 L*="A02" 

4790 GO TO 4890 
4800 Lt="AQ3" 

4810 GO TO 4890 
4820 L*="A04" 

4830 GO TO 4890 
4840 L$= "A05" 

4850 GO TO 4890 
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4860 L$= "A06 “ 

4870 GO TO 4890 
4880 L $ = " A 0 7 " 

4890 PRINT @6;L* 

4900 INPUT 86:L* 

4910 St=INT(VAL(L*)/l000) 

4920 S2=INT< <VAL(L*)-S1 *1 000 )/l00) 

4930 S3=INT<(VAL(L*)-S1*1Q00-S2*l00VI0 > 

4940 S4=VAL(L*)-S1*10Q0-S2*100-S3*l0 

4950 RETURN 

4960 PRINT "LENTER:" 

4970 FOR N=1 TO 7 

4980 GOSUB 4750 

4990 IF S4=0 THEN 5020 

5000 IF S4=7 THEN 5040 

5010 GO TO 5060 

5020 GOSUB 2310 

5030 GO TO 5200 

5040 GOSUB 2330 

5050 GO TO 5200 

5060 PRI "JCALIBRATION FACTOR <INPUT-UNITS/1 VOLT> FOR CHANNEL m jMj": -j 
5070 INPUT L 

5080 GO TO S4 OF 5090,5110,5130,5150,5170,5190 

5090 M9(M,2)=L/512 

5100 GO TO 5200 

5110 M9(M,2)=L*0.2/512 

5120 GO TO 5200 

5130 M9 < M,2)=L*5/512 

5140 GO TO 5200 

5150 M9(M,2)=L*0.1/512 

5160 GO TO 5200 

5170 M9 < M,2)=L*2/512 

5180 GO TO 5200 

5190 M9(M,2)=L*0.5/512 

5200 M9(M,1)=M 

5210 NEXT M 

5220 RETURN 

5230 PRINT UPPER CALIBRATION LEVEL FIRST!' TRY IT AGAIN!!G” 

5240 GO TO 4670 

5250 REM *********************** WRITE RECORD ************************** 
5260 DIM A(NO) 

5270 CHARSIZE 3 

5280 PRINT "LENTER FILE NUMBER (DEVICE NR.3) TO WRITE RECORD ONG : "; 

5290 INPUT D3 
5300 FIND @3:D3 
5310 PRINT @6:"P ”;D* 

5320 INPUT @6;A 
5330 A=°A-A2 (M , 2 > 

5340 WRITE 83:W*,R2,H,T1,M,M9(M,2),A2(M,2>,A 
5350 PRINT "JRECORD IS STORED IN FILE *G”;D3 
5360 DELETE A 
5370 RETURN 


6000 REM SPRGM/PART II FILE 3 

6010 REM *********************** LIN. SPECTRA ************************** 
6020 DELETE K 
6030 DIM K <CO > 

6040 A$ ="MAXIMAX SHOCK SPECTRUM” 

6050 GO TO R1 OF 6060,6080,6100,6120,6150 
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6060 K=ABS(V5> 

6070 GO TO 6170 
6080 K*ABS(G5> 

6090 GO TO 6170 
6100 K=ABS(X5) 

6110 GO TO 6170 

6120 A4="RESIDUAL SHOCK SPECTRUM” 

6130 K=ABS<Y5> 

6140 GO TO 6170 
6150 A4*"PHASE" 

6160 K=Z5 

6170 CALL "MAX”,K,Ml,14 
6180 CHARSIZE 3 

6190 PRINT "LJFREQUENCY RANGEG: 

6200 INPUT F6,F7 
6210 IF D=4 THEN 6240 

6220 VIEWPORT 8*1.79,130-4*1.79,3*2.82,100-3*2.82 
6230 GO TO 626G 

6240 VIEWPORT 8*1.79,150-4*1,79,3*2.82,100-3*2,82 
6250 PRINT @0,17:1.2,1.8 
6260 CALL "TIME”,T* 

6270 CHARSIZE 1 
6280 PRINT @D:"L";R* 

6290 PRINT @D:T* 

6300 IF R1=5 THEN 6510 

6310 IF ABS<Ml)<10 THEN 6350 

6320 IF ABS<Ml)<100 THEN 6390 

6330 IF ABS<M1)<1000 THEN 6430 

6340 IF ABS(Ml>>1000 THEN 6470 

6350 M1=*INT(M1 > + l 

6360 WINDOW F6,F7,0,M1 

6370 AXIS @D:(F7-F6)/I0,0,5 

6380 GO TO 6580 

6390 M1=INT(M1/10)*10+10 

6400 WINDOW F6,F7,0,Ml 

6410 AXIS @D:(F7-F6)/I0,5 

6420 GO TO 6580 

6430 M1=INT(M1/100)*100+100 

6440 WINDOW F6,F7,0,M1 

6450 AXIS @D:(F7-F6)/10,50 

6460 GO TO 6580 

6470 M1=INT<M1/1000)*1000+1000 

6480 WINDOW F6,F7,0,M1 

6490 AXIS @D:(F7-F6 >/I0,500 

6500 GO TO 6580 

6510 WINDOW F6,F7,0,180 

6520 AXIS @D:(F7-F6)/I0,45 

6530 MOVE @D:F6,180 

6540 DRAW @D:F7,180 

6550 DRAW 9D:F?,0 

6560 MOVE @D:(F7+F6)/2,180 

6570 GO TO 6610 

6580 MOVE 8D:F6,Ml 

6590 DRAW 9D:F7,M1 

6600 DRAW 8D:F7,0 

6610 PRINT 80,17:2,3 

6620 CHARSIZE 4 

6630 GOSUB 6700 

6640 IF Wt>0 THEN 6800 

6650 A*-W* 
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6660 PRINT 80,17:1,2,1.8 

6670 CHARSIZE 2 

6680 GOSUB 6700 

6690 GO TO 6800 

6700 IF R1=5 THEN 6730 

6710 HOVE 8D:<F7+F6>/2,Ml 

6720 GO TO 6740 

6730 MOVt: @D: (F7*F6>/2,180 

6740 PRINT 8D:"KK" j 

6750 FOR 1=1 TO LEN<A*>/2 

6760 PRINT 8D:"H"; 

6770 NEXT I 
6780 PRINT 8D:A$; 

6790 RETURN 

6800 PRINT 8D,17:1.2,2 
6810 CHARSIZE 2 

6820 H*="Natural Frequency Hz" 

6830 MOVE 8D:(F7+F6>/2,0 
6840 PRINT 8D: "JJJ"j 
6850 FOR 1=1 TO LEN<H*)/2 
6860 PRINT 8D:"H"; 

6870 NEXT I 
6880 PRINT 8D:H*; 

6890 GO TO R1 OF 6900,6940,6960,6920,6980 
6900 V$="Pseudovelocity n/sec“ 

6910 GO TO 7010 

6920 V$= "Rel. Velocity n/sec" 

6930 GO TO 7010 

6940 V*="Abs. Acceleration g" 

6950 GO TO 7010 

6960 V*= "Rel. Displacenent mm" 

6970 GO TO 7010 
6980 V*=" Phase deg" 

6990 MOVE @D:F6,90 

7000 GO TO 7020 

7010 MOVE 8D:F6,0.5*INT(Ml ) 

7020 PRINT 80:"HHHHHH"; 

7030 FOR 1=1 TO LEN<VS)/2 
7040 PRINT 80: "K " 

7050 NEXT I 

7060 FOR 1=1 TO LEN(V«) 

7070 X*=SEG(V*,I,1 ) 

7080 PRINT 80:X*;”HJ" ; 

7090 NEXT I 

7100 FOR I=F6 TO F7 STEP <F7-F6)/10 

7110 MOVE 80:1,0 

7120 PRINT 80:"HJ”;INT(I); 

7130 NEXT I 

7140 IF R1=5 THEN 7270 
7150 IF ABS(Ml)< = 10 THEN 7190 
7160 IF ABS(Ml><=100 THEN 7210 
7170 IF ABS(Ml)<=1000 THEN 7230 
7180 IF ABS<M1)>1000 THEN 7250 
7190 FOR 1=0 TO ABS(Ml) STEP 1 
7200 GO TO 7280 

7210 FOR 1=0 TO ABS(Ml) STEP 10 
7220 GO TO 7280 

7230 FOR 1=0 TO ABS(Ml) STEP 100 
7240 GO TO 7280 

7250 FOR 1=0 TO ABS(Ml) STEP 1000 








7260 GO TO 7280 

7270 FOR 1=0 TO 180 STEP 45 

7280 HOVE 8D:F6,I 

7290 PRINT 8D:"HHH" jI; 

7300 NEXT I 
7310 MOVE BD:F6,0 
7320 PRINT 80,17:1 ,2 
7330 A$="K(Zero Danping)" 

7340 PRINT 80:"JJJJ";A« 

7350 1=1 

7360 FOR 1=1 TO CO 

7370 IF F5(I)> 0.01 THEN 7390 

7380 NEXT I 

7390 MOVE 8D:F5(I>,ABS(K(I)> 

7400 FOR 1=1 TO CO 

7410 IF F5(I>=0.01 THEN 7360 

7420 DRAW 80:F5(I),ABS(K<I) ) 

7430 NEXT I 
7440 RETURN 

7450 REM ***************** WRITE DATA ft********************************: 
7460 CHARSIZE 3 

7470 PRINT "JLGENTER FILE NUMBER (DEVICE NR.3) TO WRITE DATA ON: "; 

7480 INPUT 01 
7490 FIND S3;D1 

7500 WRITE 83:CO,W*,F5,X5,V5,G5,T5,Y5 > Z5 
7510 PRINT "RESULTS STORED IN FILEG *" jDl 
7520 RETURN 

7530 REM ****************** READ DATA ********************************* 
7540 CHARSIZE 3 

7550 PRINT "JLGENTER FILE NUMBER (DEVICE NR.3) TO READ DATA FROM: "j 
7560 INPUT D3 

7570 DELETE F5,X5,V5,G5,T5,Y5,Z5,K,A 
7580 FIND 83:D3 
7590 READ e3:C0 
7600 FIND 83:D3 

7610 DIM F5(C0),X5(CQ),V5(C0),G5(C0),T5(C0),Y5(C0),Z5(CO),K(CO) 

7620 READ 83:CO,W*,F5,X5,V5,G5,T5,Y5,Z5 

7630 IF W1>0 THEN 7650 

7640 PRINT "DATA ARE IN MEMORY!G* 

7650 RETURN 

7660 REM *********************** LOG. SPECTRA ********************** 
7670 VIEWPORT 5*1.79,110-4*1.79,10+3*2.82,100-3*2.82 
7680 GO TO 7740 

7690 VIEWPORT 8*1.79,130-4*1.79,3*2.82,100-3*2.82 
7700 IF D=4 THEN 7720 
7710 GO TO 7740 

7720 VIEWPORT 8*1.79,150-4*1.79,3*2.82,100-3*2,82 

7730 PRINT 8D,17:1.2,1.8 

7740 CHARSIZE 1 

7750 DELETE K 

7760 DIM K(CO) 

7770 CALL "TIME",T$ 

7780 PRINT 8D:"L " ;R* 

7790 PRINT 8D:Ti; 

7800 REM 

7810 WINDOW LGT(N6),LGT(N7),LGT(N1),LGT<N2) 

7820 IF P >10 THEN 7980 
7830 HI=INT(LGT(N1>) 

7840 FOR 1 = 10 A H1 TO 10MHI+1) STEP 10"H1»P 
7850 MOVE 8D:LGT(N6),LGT(I) 











7860 DRAW BD:LGT(N7),LGT(I) 

7870 NEXT I 
7880 HI=H1>1 

7890 IF HI<=LGT<N2> THEN 7840 
7900 HI=LGT<N6) 

7910 FOR 1 = 10 A H1 TO 10 A <H1+1> STEP 10 A H1*P 
7920 MOVE SD:LGT< I) ,LGT(N1> 

7930 DRAW SD:LGT<I),LGT(N2) 

7940 NEXT I 
7950 H1=H1+1 

7960 IF HI<=LGT<N7) THEN 7910 

7970 GO TO 8020 

7980 AXIS SD;1,1,0,-2 

7990 MOVE BD:LGT<N6),LGT(N2) 

8000 DRAW BD:LGT(N7),LGT(N2) 

8010 DRAW BD:LGT<N7),LGT<N1> 

8020 PRINT BD,17:2,3 
8030 CHARSIZE 4 

8040 GO TO R OF 8050,8070,8090 
8050 A$="MAXIMAX SHOCK SPECTRUM" 

8060 GO TO 8100 

8070 A*="RESIDUAL SHOCK SPECTRUM" 

8080 GO TO 8100 

8090 A$="MAXIMAX and RESIDUAL SHOCK SPECTRUM" 

8100 GOSUB 8170 

8110 IF W1>0 THEN 8240 

8120 A$=W* 

8130 PRINT BD,17:1.2,1.8 
8140 CHARSIZE 2 
8150 GOSUB 8170 
8160 GO TO 8240 

8170 MOVE BD:LGT(10 A < 0.5*LGT < N6»N7))),LGT(N2 > 
8180 PRINT BD:"KK“ ; 

8190 FOR 1=1 TO LEN(A«)/2 
8200 PRINT BD:"H"; 

8210 NEXT I 
8220 PRINT BD:A*; 

8230 RETURN 

8240 GO TO R OF 8250,8270,8290 
8250 V$="Pseudovelocity n/sec" 

8260 GO TO 8300 

8270 V*="Rel, Velocity «/sec" 

8280 GO TO 8300 

8290 V«=”Velocity n/sec " 

8300 MOVE BD:LGT(N6),LGT(10 A < 0.5#LGT(N2*N1))) 

8310 PRINT BD,17:1.5,2.5 

8320 CHARSIZE 2 

8330 PRINT BD:“HHHHHHH"j 

8340 FOR 1*1 TO LEN<vi>/2 

8350 PRINT BD: ”K"; 

8360 NEXT I 

8370 FOR 1=1 TO LEN(V») 

8380 X*=SEG<V*,I,1) 

8390 PRINT 8D:X*;"HJ"; 

8400 NEXT I 

8410 MOVE BD:LGT(10 A < 0.5»LGT < N6*N7))>,LGT(N1> 
8420 H$="Natural Frequency Hz" 

8430 FOR 1*1 TO LEN<H«)/2 
8440 PRINT BD:"H" ; 

8450 NEXT I 
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8460 PRINT 8D:“JJJ", H$; 

8470 HI=0 

8480 FOR I=LGT(N1> TO LGT<N2> STEP 10 A H1 
8490 MOVE @D:LGT<N6),I 

8500 PRINT 8D;"HHHHH"; " 1 0 A " ;INT< HI+LGT<N1)) 
8510 H1=H1+1 
8520 NEXT I 
8530 HI=0 

8540 FOR I=LGT(N6) TO LGT < N7) STEP 1Q A H1 
8550 MOVE 8D: I,LGT<N1) 

8560 PRINT @D:"HHJ";"10 A “;INT(H1+LGT<N6>) 
8570 H1=H1+1 
8580 NEXT I 

8590 MOVE @D:LGT < N6) ,LGT < N1) 

8600 IF D=32 THEN 8620 
8610 PRINT @0,17:1,2 
8620 CHARSIZE 2 
8630 A*="(Zero Dinping)" 

8640 PRINT 8D:”JJJ";A* 

8650 GOSUB 8770 
8660 GOSUB 8820 
8670 IF R=3 THEN 8690 
8680 RETURN 
8690 R=R+1 

8700 IF 004 THEN 8740 

8710 PRINT "LYOU MAY CHANGE THE PEN COLOR G* 
8720 PRINT ■ Press RETURN to continue" 

8730 INPUT A* 

8740 GOSUB 8770 
8750 GOSUB 8820 
8760 RETURN 

8770 GO TO R OF 8780,8800,8780,8800 
8780 K=V5 
8790 GO TO 8810 
8800 K=Y5 
8810 RETURN 
8820 1=1 

8830 FOR 1=1 TO CO 
8840 IF F5<I)>0.01 THEN 8870 
8850 NEXT I 
8860 RETURN 

8870 MOVE @D:LGT < F5(I)),LGT < ABS < K<I) ) ) 

8880 FOR 1=1 TO CO 
8890 IF F5<I)=0,01 THEN 8830 
8900 IF R<*3 THEN 8920 
8910 DASH 170 

8920 DRAU 80: LGT < FS<I )) ,LGT<ABS (K<I) ) ) 

8930 NEXT I 

8940 IF W1>0 THEN 8960 
8950 DASH 0 
8960 HOME 
8970 RETURN 

8980 REM »*»#**»**»**»« NUMERICAL OUTPUT 
8990 WINDOW 0,130,0,100 
9000 VIEWPORT 0,130,0,100 
9010 11=0 
9020 12=1 
9030 N-l 

9040 IF D=32 THEN 9070 
9050 IF D=40 THEN 9130 












9060 PRINT 80,17:1.5,2.5 
9070 CHARSIZE 3 
9080 PRINT 80:'%" 

9090 HOVE 80:65,99 
9100 FOR 1=1 TO LEN(U*>/2 
9110 PRINT 80:"H"; 

9120 NEXT I 
9130 PRINT 8D:W* 

9140 CHARSIZE 1 
9150 CALL “TIME",T$ 

9160 PRINT 80: "*J“ ;R*, “ ,, jT*,"JI H 
9170 IF 0=40 THEN 9190 

9180 PRINT 80: USING 9190:"KKK"PAGE12;"JJ“ 

9190 IMAGE 3A,122X,4A,20,2A 
9200 CHARSIZE 2 

9210 PRINT 80:" FREQ SHOCK SPECTRUM FOURIER TRANS" 

9220 PRINT 80:" <HZ) <MM) <M/S> (G) <MS> <M/S) (DEG>_“ 

9230 IMAGE 40.1D,X,50.20,X,40.20,X,50.ID,X,40.20,X,40.20,2X,40 
9240 IMAGE 65T,40.1D,X,50.20,X,40.20,X,50,1D,X,4D.20,X,40.20,2X,4D 
9250 FOR 1=1+11 TO CO 
9260 IF I>50»N THEN 9300 

9270 PRINT 80: USING 9230:F5<I>,X5<I),V5<I>,G5(I>,T5(I),Y5(I>,Z5<I) 

9280 NEXT I 
9290 RETURN 

9300 IF 0=40 THEN 9350 

9310 CHARSIZE 1 

9320 PRINT eD:“ A JJJJ2" 

9330 CHARSIZE 2 
9340 PRINT 80:"J" 

9350 FOR I=1+50*N TO CO 
9360 IF I>(N+l)*50 THEN 9430 
9370 IF 0=40 THEN 9400 

9380 PRINT 80: USING 9240:F5<I),X5(I),V5(I>,G5<I),T5(I>,Y5(I),Z5(I) 

9390 GO TO 9410 

9400 PRINT 80: USING 9230:F5<I),X5<I),V5<I>,G5<I),T5CI>,Y5<I),Z5(I> 

9410 NEXT I 

9420 RETURN 

9430 IF 0=40 THEN 9450 

9440 CALL "UAIT",4 

9450 11*100 

9460 12*12+1 

9470 N-N+2 

9480 IF 0=40 THEN 9500 
9490 GO TO 9040 
9500 GO TO 9250 

9510 REM «•»*»»*»•*»*•***•»**» ANALYSIS ft******************************* 

9520 R»3 

9530 GOSUB 7670 

9540 R*1 

9550 MOVE 80:3.11,2 

9560 PRINT 80: USING "2A,5X,4A,3X,7A":"NR","CHZ[M/SEC]" 

9570 POINTER X,Y,Q* 

9580 IF Q*«"0" THEN 9750 

9590 MOVE 80:X-LGT(1,022>,Y+LGT(1.022) 

9600 CHARSIZE 2 

9610 PRINT 80:“J*" jR 

9620 GOSUB 9650 

9630 R-R+l 

9640 GO TO 9570 

9650 VIEWPORT 106,130,0,100 
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9660 WINDOW 0,32,0,100 
9670 HOVE 80:0,94 
9680 FOR 1*1 TO R 
9690 PRINT 8D:"JJ" \ 

9700 NEXT I 

9710 PRINT 8D: USING 9720:”J*” ,R,“:“,13*X,10 A Y 
9720 IMAGE 2A,20,A,2X, 4D .1D,3X,3D.D 

9730 VIEWPORT 5*1.79,110-4*1.79,10+3*2.82,100-3*2.82 
9740 WINDOW LGT(N6),LGT(N7),LGT<N1),LGT<N2) 

9750 DASH 0 
9760 HOME 
9770 RETURN 

9780 REM ********************* MULTI DISPLAY *************************** 
9790 CHARSIZE 3 
9800 D-32 

9810 PRINT “LJTENTER NUMBER OF DESIRED SPECTRAG: “; 

9820 INPUT wT 
9830 DIM W2(W1) 

9840 PRINT " ENTER FILE NUMBERS ON DEVICE NR.3G: "; 

9850 INPUT W2 

9860 PRINT “_ENTER: 1 FOR MAXIMAX-SPECTRA_2 FOR RESIDUAL-SPECTRAg: "> 

9870 INPUT R 
9880 D3*W2<1> 

9890 GOSUB 7570 

9900 GOSUB 7690 

9910 COSUB 10060 

9920 IF Wl-1 THEN 10050 

9930 FOR 10-2 TO W1 

9940 D3=W2(I0) 

9950 GOSUB 7570 
9960 DASH 11*10 
9970 GOSUB 8770 
9980 GOSUB 8820 
9990 GOSUB 10060 
10000 NEXT 10 
10010 DASH 0 
10020 Wl-0 

10 030 DELETE W2,X,Y,10,D3,Q* 

10040 W*-"TEST" 

10050 RETURN 

10060 POINTER X,Y,Q* 

10070 MOVE 8D:X,Y 
10080 DRAW 8D:X+Q,25,Y 
10090 PRINT BD;" "jW* 

10100 RETURN 
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1 REM ///////////////// 8SPRGM/PART III \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\ 

2 INIT 

3 GO TO 260 

4 IF 0*3 THEN 94 

5 GO TO 7450 

8 IF 0=3 THEN 94 

9 GO TO 7530 

12 IF 0=3 THEN 94 

13 GO TO 8980 
16 GO TO 190 
20 GO TO 100 

24 IF 0=3 THEN 94 

25 R = 1 

26 GO TO 7690 

28 IF 0=3 THEN 94 

29 R*2 

30 GO TO 7690 

32 IF 0=3 THEN 94 

33 R=3 

34 GO TO 7690 

36 IF 0=3 THEN 94 

37 GO TO 9510 

40 IF 0=3 THEN 94 

41 CO TO 9780 
45 RETURN 

48 IF 0=2 THEN 90 

49 GO TO 510 

50 DELETE A 

52 IF 0=2 THEN 90 

53 GO TO 1780 

54 DELETE A 

64 IF 0=3 THEN 94 

65 Rl=l 

66 GO TO 6010 

68 IF 0=3 THEN 94 

69 R1=4 

70 GO TO 6010 

72 IF 0=3 THEN 94 

73 R1*2 

74 GO TO 6010 

76 IF 0=3 THEN 94 

77 R1=3 

78 GO TO 6010 

80 IF 0=3 THEN 94 

81 R1=5 

82 GO TO 6010 

90 PRINT "You got to PRESS KEY 5 first'C" 

92 RETURN 

94 PRINT "You got to PRESS KEY 4 firstIG" 

96 RETURN 

100 REM»*«**»*»******** LOAD ROUTINE a********************************** 
110 IF 0=3 THEN 170 

120 DELETE F5,X5,V5,G5,T5,Y5,Z5,C,C0,K,A 
130 C0=90 
140 0=3 

150 APPEND "8SPRGM/PARTIII"j1 

160 PRINT "JGSPRGM/PART III IS LOADED!" 

170 DELETE 5140,20000 
180 RETURN 
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190 IF 0-2 THEN 240 
200 0-2 

210 DELETE 270,20000 

220 APPEND "8SPRGM/PARTII";260 

230 RENUMBER 6000,10,260 

240 PRINT "JGSPRGM/PART II IS LOADED!" 

2S0 RETURN 

260 REM************** DEFAULT CONDITION ******************************** 
270 W*= "TEST“ 

280 0-3 
290 Ul-0 
300 D-32 
310 Fl-1 
320 F2-1000 
330 F8-1000 
340 Nl-0.01 
350 N2-100 
360 N6-1 
370 N7-1000 
380 CO-90 
390 N0-4096 
400 Kl-2 
410 K 2-2 
420 K3-10 
430 P-10 

440 R*-“NRL - CODE 5842" 

450 DIM M9(7,2>,A2(7,2),C(C0> 

460 M9=0 
470 A2-0 

480 DIM D»<8>,H*(20>,V*(22>,G*<7>,M*<7>,N«(8>,X*(22>,PS(35>,L*(4>,Q*(1 > 
490 PRINT “ GSPRGM/PART III IS LOADED (DEFAULT CONDITIONS)!" 

500 RETURN 

510 REM A**************** READ RECORD*********************************** 
520 CHARSIZE 3 

530 PRINT "LJENTER FILE NUMBER (ON DEVICE 3) TO READ RECORD FROMG: "; 

540 INPUT D1 
550 DIM A(NO) 

560 U3-130 

570 IF D< >4 THEN 600 
580 U3-150 

590 PRINT BD,17:1.2,1.8 
600 CHARSIZE 2 
610 CALL "TIME", TO 
620 PRINT BD: "l*K",R* 

630 PRINT SD:T* 

640 FIND S3:D1 

650 READ B3:U«,R2,H,T1,M 

660 FIND B3:D1 

670 READ 83:W*,R2,H,T1,M,M9(M,2>,A2(M,2>,A 
680 M9(M , 1)-M 
690 A2(M,1)-M 
700 S-M9(M,2> 

710 WINDOW 0,W3,0,100 
720 VIEWPORT 0,W3,0,100 
730 IF W1>0 THEN 810 
740 PRINT BD,17:2,3 
750 CHARSIZE 4 
760 MOVE BD;W3/2,98 
770 FOR 1-1 TO LEN(W*>/2 
780 PRINT BD :"H”; 
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790 NEXT I 
800 PRINT @D:U« 

810 PRINT @0,17:1.2,1.8 
820 CHARSIZE 2 

830 GO TO R2 OF 840,970,1110 

840 GOSUB 1280 

850 S1=S 

860 GOSUB 1460 

870 PRINT @D:“*" 

880 CALL “DIF3",A,A 

890 S1=S1/H 

900 GOSUB 1340 

910 GOSUB 1460 

920 CALL "DIFS",A,A 

930 S1=S1*10 A 3/(H*9.8044) 

940 GOSUB 1400 
950 GOSUB 1460 
960 GO TO 1230 
970 GOSUB 1340 
980 S1=S 
990 COSUB 1460 
1000 PRINT @D: 

1010 CALL "INT",A,A 

1020 S1=S1*H 

1030 GOSUB 1280 

1040 GOSUB 1460 

1050 CALL "DIF3”,A , A 

1060 CALL "DIF3",A,A 

1070 Sl a Sl*l0 A 3/(H A 2#9.8044) 

1080 GOSUB 1400 
1090 GOSUB 1460 
1100 GO TO 1230 
1110 GOSUB 1400 
1120 S1=»S 
1130 GOSUB 1460 
1140 PRINT @D: 

1150 CALL "INT",A,A 

1160 S1=S1*1.0E~3*H*9.8044 

1170 GOSUB 1340 

1180 GOSUB 1460 

1190 CALL "INT",A,A 

1200 S1=S1»H 

1210 GOSUB 1280 

1220 COSUB 1460 

1230 VIEWPORT 0,W3,0,100 

1240 WINDOW 0,W3,0,100 

1250 HOVE @D:0,0 

1260 PRINT @D:"* orig. tiMe record"; 
1270 RETURN 

1280 VIEWPORT 5,W3-30,4,32 
1290 V*»"DISPLACEMENT" 

1300 G*-"D :" 

1310 M*»"D ciin : " 

1320 N9»" C mm 1" 

1330 RETURN 

1340 VIEWPORT 5,W3-30,36,64 
1350 V**"VELOCITY" 

1360 G**“V mix :" 

1370 M*-"V Min :" 

1380 N««" Cm/scc1 * 










1390 RETURN 

1400 VIEWPORT 5,W3-30,68,96 
1410 VS*"ACCELERATION" 

1420 GS*"A Max :" 

1430 MS*"A Min 
1440 NS*" tgl" 

1450 RETURN 

1460 CALL "MAX",A,VI,14 
1470 CALL "MIN",A,V2,13 
1480 WINDOW 0,N0,V2,V1 
1490 IF D*32 THEN 1550 
1500 MOVE SD:0,0 
1510 FOR 1*1 TO NO 
1520 DRAW 8D:I,A(I) 

1530 NEXT I 

1540 GO TO 1560 

1550 CALL "DISP",A 

1560 AXIS 8D:400,0 

1570 FOR 1*0 TO NO STEP 2000 

1580 MOVE @D:I,0 

1590 PRINT 8D:"HJ“,I*H; 

1600 NEXT I 

1610 PRINT @D:"HHHHHHHHHHHIMS]" 

1620 MOVE @D:420 0,vT 

1630 PRINT SDi USING 1670:"J",GS,V1*S1,NS; 

1640 MOVE eD:4200,VI 
1650 PRINT eD!"J "} 

1660 PRINT 8Di USING 1670!"J",MS,V2*S1,NS; 

1670 IMAGE A,7A,6D.1D, 11A 
1680 MOVE 8D:0,0.5*<V2+V1) 

1690 PRINT eD:”HH"; 

1700 FOR 1*1 TO LEN(VS)/2 
1710 PRINT 8D:"K "; 

1720 NEXT I 

1730 FOR 1*1 TO LEN<VS) 

1740 XS*SEG<VS,1,1) 

1750 PRINT eD:XS;"HJ"j 
1760 NEXT I 
1770 RETURN 

1780 REM ******************** PROCESS X********************************* 
1790 CHARSIZE 3 
1800 CALL "TIME",BS 
1810 PRINT 8Di "l-'jBSl " 1 " 

1820 PRINT 832,26:2 
1830 IF FI<1 THEN 1890 
1840 IF FI>1000 THEN 1890 
1850 IF F2<1 THEN 1900 
1860 IF F2>1000 THEN 1900 
1870 IF F2<*F1 THEN 1900 
1880 GO TO 1930 
1890 IF C0=1 THEN 1930 

1900 PRI "ANALYSIS FREQUENCY IS OUT OF RANGE ENTER NEW VALUES OF F1,F2G" 
1910 INPUT FI,F2 
1920 GO TO 1800 

1930 DELETE F5,XS,V5,G5,T5,Y5,Z5,C,K , A 
1940 DIM A<NO) 

1950 C3*INT<(MEMORY/8-18>/8> 

1960 DIM F5(C0>,X5(C0>,V5(C0>,G5(C0>,T5(C0>,Y5(C0>,Z5(C0> ,C(C0> 

1970 FIND S3:D1 

1980 READ S3:WS,R2,H,T1,M,S 
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1990 FIND @3:01 

2000 E1=INT(1000/<K1*F8*H>> 

2010 £2 = INT < 1000/<K2*F0»H)) 

2020 E3=INT(1000/(K3*F8*H)) 

2030 F5(1)=F1 

2040 IF C0 = 1 THEN 2080 

2050 FOR 1=2 TO CO 

2060 F5<I)=<F2/F1) A <1/<C0-1)>*F5<1-1> 

2070 NEXT I 
2080 X5=0.01 
2090 V5=0.01 
2100 G5=0.01 
2110 T5=0.01 
2120 Y5= 0.01 
2130 Z5=0.01 
2140 11=1 
2150 12=0 

2160 IF FI<10 THEN 2200 
2170 12=1 

2180 IF FI<100 THEN 2200 

2190 12=2 

2200 GO TO 2370 

2210 E=E1 

2220 GOSUB 2460 

2230 IF I1>C0 THEN 2310 

2240 GO TO 2370 

2250 E=E2 

2260 GOSUB 2460 

2270 IF IOCO THEN 2310 

2280 GO TO 2370 

2290 E=E3 

2300 GOSUB 2460 

2310 DELETE A 

2320 CALL "TIME",T* 

2330 GOSUB 3320 
2340 PRINT @32,26:0 

2350 PRINT @D:"_HAX. POSSIBLE ANALYSIS FREQUENCIES 
2360 RETURN 
2370 F3=F5<11) 

2380 C1=I1 

2390 12=12+1 

2400 FOR 11=1 TO CO 

2410 IF F5(I1)>10 A I2 THEN 2430 

2420 NEXT II 

2430 F4=F5<11-1 ) 

2440 GO TO 12 OF 2210,2250,2290 
2450 RETURN 

2460 READ @3:U*,R2,H,T1,M,S,A2<M,2>,A 

2470 FIND @3:D1 

2480 IF E< =1 THEN 2550 

2490 J=1 

2500 FOR 1=1 TO NO STEP E 
2510 A(J)=A(I> 

2520 J=J+1 
2530 NEXT I 
2540 GO TO 2560 
2550 E=1 
2560 N5=N0/E 
2570 H2=H*E 

2580 GO TO R2 OF 2590,2630,2660 







2590 CALL "DIF3",A ( A 

2600 CALL “DIF3", A,A 

2610 S2=S 

2620 GO TO 2670 

2630 CALL ”DIF3",A,A 

2640 S2=S*H2 

2650 GO TO 2670 

2660 S2=»S*10*-3*H2*2*9,8044 

2670 C2=I1-1 

2680 P1=0.0062831853*F3*H2 
2690 FOR I=C1 TO C2 
2700 C(I)=2*C0S(PI) 

2710 IF Cl=C2 THEN 2730 
2720 P 1 =P 1 •*< F4/F3) A (1/<C2-C1>) 

2730 NEXT I 

2740 IF P1< 0.902 THEN 2790 
2750 IF 143.566/H2>F2 THEN 2790 

2760 PRINT "LGUARNINGi FREQUENCIES ARE TOO HIGH FOR YOUR TIME INCREMENT" 
2770 PRINT " ERRORS ARE OUTSIDE THE RANGE -7 TO +3 PERCENT FOR" 

2780 PRINT " FREQUENCIES ABOVE '*; 143.566/H2 ; ” HERTZ” 

2790 FOR I=C1 TO C2 
2800 F=F5(I) 

2810 IF F< = 1 THEN 2840 
2820 N=1000/(H2*5*LGT<F>) 

2830 IF N<N5 THEN 2850 

2840 N=N5 

2850 X1=0 

2860 X2=0 

2870 R3=0 

2880 R4=0 

2890 IF F=0 THEN 2920 
2900 N4=N-1-500/<F*H2) 

2910 IF N4 >3 THEN 2930 
2920 N4*3 

2930 FOR J*1 TO N-l 
2940 X3=C(I)*X2-X1+A<J) 

2950 IF J<N4 THEN 3000 

2960 IF ABS(X3)<A6S(R4) THEN 3030 

2970 R4»X3 

2980 T4*J 

2990 GO TO 3030 

3000 IF ABS(X3)<ABS<R3) THEN 3030 

3010 R3=X3 

3020 T3=J 

3030 XI=X2 

3040 X2*X3 

3050 NEXT J 

3060 IF ABS(R3)>ABS(R4) THEN 3090 

3070 R3=R4 

3080 T3=T4 

3090 X=S2*R3 

3100 V=0.00628318*F*X 

3110 G=6.288318*F#V/9,8044 

3120 T=Tl+(T3-0.5>*H2 

3130 Y=0.00628318*F»S2*R4 

3140 Z=l.0E-3*(T1+(T4-0.5)*H2)*F+0.25 

3150 Z=Z-INT(Z) 

3160 IF Y >0 THEN 3180 
3170 Z=Z-0.5 

3180 IF Z< *0,5 THEN 3200 
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3190 Z=Z-1 
3200 Z=360*Z 
3210 F5(I)=F 
3220 X5(I)=X 
3230 U5(I)=V 
3240 G5(I)»G 
3250 T5(I>=T 
3260 Y5(I)=Y 
3270 Z5(I)=Z 

3280 PRINT @D: USING 3290;I,F,U,N; 

3290 IMAGE 20T,3D,4X,4D.1D,4X,5D,2D,3X,6D 
3300 NEXT I 
3310 RETURN 

3320 REM **#*»•* SUBROUTINE TIME ***********»*»***#*###***##*##»x##***# 
3330 DIM B(2) 

3340 C$=B* 

3350 FOR 1=1 TO 2 
3360 X*=SEG(C*,U ,8) 

3370 Y*=SEG(X*,4,2> 

3380 Z4=SEG(X*,7,2) 

3390 B(I)=VAL(X*> + (VAL(Y*>*60+VAL(Z*> >/360 0 
3400 C*=T* 

3410 NEXT I 
3420 U=B(2)-B(1) 

3430 B3=(U-INT(U))*60 
3440 B4=(B3-INT(B3))*60 

3450 PRINT @D:“..CALCULATION TIME <hf:win:sec); 

3460 PRINT @D: USING 3470:INT<U),":“,INT(B3),“:",INT(B4) 

3470 IMAGE 2D,A,2D,A,2D 
3480 RETURN 

5000 REM - 

5010 REM - 

5020 REM THIS PROCEDURE MAY BE USED TO CALCULATE SHOCK SPECTRA OF 

5030 REM TIME RECORDS STORED ON EXT. TAPE (FILE 1 TO 7) AND TO STORE 

5040 REM THE ACCORDING DATA ON INT. TAPE (FILE * TO FILE *+7>. 

5050 FOR Dl=l TO 7 
5060 GOSUB 52 
5070 N=D1+0 
5080 FIND N 

5090 MARK 1,7000 

5100 FIND N 

5110 WRITE C0,W*,F5,X5,V5,G5,T5,Y5,Z5 
5120 NEXT D1 
5130 RETURN 
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